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We report about the current status of our ongoing study of the chiral limit of two-flavor QCD at 
finite temperature with twisted mass quarks. We estimate the pseudo-critical temperature Tc for 
three values of the pion mass in the range of mps c^; 300 and 500 MeV and discuss different chiral 
scenarios. 

Furthermore, we present first preliminary results for the trace anomaly, pressure and energy den- 
sity. We have studied several discretizations of Euclidean time up to A^^ = 12 in order to assess 
the continuum limit of the trace anomaly. From its interpolation we evaluate the pressure and 
energy density employing the integral method. Here, we have focussed on two pion masses with 
mps ~400and 700 MeV. 
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1. Introduction 

The order of the phase transition in the case of two-flavor QCD in the chiral limit remains 
an open question. While with Wilson quarks the transition is found to be compatible with 

a second order phase transition in the universality class of a 3d 0(4) spin model, a first order 
transition seems to be favoured in analyses with staggered fermions at A^^ = 4 [|3|, 

The thermal equation of state (EoS) constitutes a relevant ingredient in the hydrodynamic 
evolution of the quai^k-gluon plasma created in heavy-ion experiments. It can be determined non- 
perturbatively in lattice calculations. In the recent past the EoS has been studied extensively using 
the staggered type of quark discretization, mostly with Nf = 2 + \ flavors at the physical point 
[Q, ^. The much more compute-intensive Wilson-like discretizations are less investigated however 
[^, ^]. In the latter study the fixed scale approach is used as compared to the more traditional fixed 
Nx approach. 



2. Lattice Setup 

The lattice setup in our ongoing investigations equals the one employed by the European 
Twisted Mass Collaboration (ETMC) for their Nf = 2 simulations [^. It employs the twisted 
mass action in terms of the twisted fields % = exp {—inysTT, / 4)\if 

Sf[U, Y, V] = (4,, - KDw{x,y)[U]+2iKaHY5^35,,y)x{y) ■ (2-1) 

.v.v 

in the quark sector, while the gauge sector is described by the tree-level Symanzik improved gauge 
action 

Sfy-^[U] = p (co£[l - -ReTr (t/p)] +C1 ^[1 - -ReTr ([/«)]) . (2.2) 

The latter two sums extend over all possible plaquettes (P) and all possible planar rectangles (R), 
respectively. 



3. Pseudo-Critical Temperatures and Chiral Limit 

For the present study of the chiral Umit we rely on simulations with A^^- = 12 at pion masses 



mps ~ 320 MeV, 400 MeV and 470 MeV that have been analyzed in Ref. [|10p (for historical 
reasons we call these ensembles A12, B12 and C12). Our determination of the pseudo-critical 
temperature is based on the measurement of the variance of over the gauge ensemble 

4v = ^(((w)'>-(w)') • (3.1) 

It corresponds to the disconnected part of the usual chiral susceptibility and should show a maxi- 
mum in the region of Tc. This is indeed the case for all our ensembles and two representative cases 
are shown in the two left panels of Fig. [T[ From fitting Gaussian function to the data of a? ^ around 
the maxima we infer values of the pseudo-critical couplings pc that are converted to a physical 
value of Tc using an interpolation of a(j8) JTO[]. At leading order in chiral perturbation theory and 
for a phase transition of second order the pion mass dependence of is expected to be given as 

Tc{m^) = UO) +A inj^^^^ , (3.2) 
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where rt (0) is the critical temperature in the chiral limit and p and 5 are critical exponents cor- 
responding to the universality class of second order phase transition. We have restricted our- 
selves to the chiral scenarios discussed in Ref. [ |l0| ] including a first order scenario as well as 
the 0(4) and Z(2) second order scenarios, for the latter assuming a second order endpoint located 
at rrijc^c = MeV or alternatively at m^^.c = 200 MeV. The result of fits of Eq. (3^) to our data 
is shown in the right most panel of Fig. |[ As the fitted curves are all describing the given data 
quite well we conclude that the present set of pion mass values can not discriminate among the 
different chiral scenarios that have been studied. For the 0(4) model the fit prefers a value of 
Tc{0) = 152 (26) MeV. 
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Figure 1: Determination of Tc from using a Gaussian fit to the maximum for the A12 and B12 ensem- 
bles. Comparison of different scenarios for the chiral limit of Tc. 



4. The Trace Anomaly 

For the EoS we concentrate on one of the values of pion masses used in the chiral limit study 
above, namely the one corresponding to nips ~ 400 MeV. We have added at the same pion mass 
additional runs at smaller Nt = 4, 6 and 8 (henceforth denoted by B4, B6, B8 . . .). Apart from the 
latter, for which A^^ = 28 has been chosen, all lattices have a spatial extent No = 32. Moreover, 
further ensembles at mps ~ 700 MeV (further on referred to as the D mass) were generated with 
sizes nIxNt: = 24^ x 10, 20^ x 8 and 16^ x 6 (referred to as DIO, D8 and D6). 

from deriva- 



The direct evaluation of pressure p = T 



and energy density £ = ^ 



V 

tives of the partition function Z is problematic given the lattice spacing dependence of both the 
temperature T = 1/ (N^a) and the volume V = N^a^. The by now standard approach is the use of 
the integral method to calculate the pressure as a temperature integral of the total derivative of the 
partition function with respect to the lattice spacing, the so called trace anomaly: 

/ _ s-3p _ _ T /d\nZ\ 
~ ~ ~Vr4 \ Jlna/^^b 

= <^/^]^| j(^^Tr£[/p\ +^/ReTr£[/«\ +B, (^Dw[C/]z)sub (4.1) 

I \ ^ / sub \ ^ / sub 



sub 
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Here, Bp, and B^ are related to the jS-functions, the derivatives of the bare pai^ameters with 
respect to the lattice spacing, as follows: 

It is necessary to subtract from each term in above expression the corresponding T = vacuum con- 
tribution (. . = (. . .)7->o ~ (• • •)r=o order to achieve a finite result. Details of the subtraction 
on the basis of the available T = lattice data as well as on the evaluation of the /3 -functions will be 
given in the following section. In the left panels of Figs. |3| and ^ we show the trace anomaly for the 
two cases of pseudoscalar masses under investigation. In both cases we observe sizeable lattice ar- 
tifacts in the height of the maximum and even in the falling edge at larger temperatures. Moreover, 
the precision in case of the smaller mass is not yet satisfactory, especially at small temperatures. 

5. /3 -Functions and T = Subtraction 

For evaluating the three j3-functions we consider fits to lattice data of the Sommer scale ro in 
the chiral limit (denoted by {'-f))- To this end the correct asymptotic behavior is built into the fit 
functions explicitly following Ref. [jlT]]. For instance we determine Bp via the identity 

by fitting r^/ato the formula 

\a/ dQ{a2L{p) + diR{l5y) a2L(Psub) 

The ratio R{P) is defined in terms of the known two-loop perturbative formula a2L(j8) and jSsub = 
3.9 has been chosen in above formula. The three parameter fit of Eq. (5.2) to (— ) (see the left 




panel of Fig. ^ yields ;^^/dof = 1.2. The thus obtained j3-function is shown in the middle panel 
of Fig. |2[ The interpolation provided by the fit of Eq. ( |5.2| ) has also been used to set the scale using 



the physical value of ro = 0.420(15) fm by ETMC 



The second j8 -function associated with the mass is evaluated from a similar identity [11]: 



r^jj. as well as its derivative are obtained by fitting the following expression to r;^./^: 

r,At=H^ PiP), P{P) = a,{l+b,R{Pf) , (5.4) 



-r) 

where jSq = [11 -IN f /?>)/ [An f and /o = l/(27r2). The third and remaining j8 -function involving 
Bk is calculated in the most straight-forward manner from an explicit derivative of jq. with respect 



to j8 using the Pade interpolation of Ref. [13] 
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In order to obtain (.•.)sut>, i- e. to subtract the T = expectation values, we have used all 
available lattice data from ETMC. For these it has been necessary to interpolate in a/i using spline 
functions to match with the simulated bare mass at T > 0. Further additional T = runs have been 
simulated in order to perform the subtraction more reliably. However, not all simulation points at 
finite temperature are supplemented with an associated T = simulation. Thus, we have performed 
an interpolation in j8 using a polynomial ansatz of fifth order. For the plaquette, the rectangle and 
the Wilson hopping term Dw we have obtained values for per degree of freedom of 2.7, 2.3 
and 2.9, respectively. The remaining term, for which no fit of reasonable quality could be obtained 
using this ansatz, has been interpolated using splines. For the D ensembles with larger mass, 
sufficient T = data newly generated is available. Hence, no interpolations in the bare coupling 
and only few interpolations in the bare mass had to be done. 
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Figure 2: Left: Interpolation of ja in the bare coupling. The point at j3 = 4.6 is not obtained from a chiral 



extrapolation and is not included in the fit. Middle: The j3-function obtained according to Eq. (5.1 ). We also 
show the perturbative 2-loop expectation at large couplings (j3) = — 12j3o — 72-^. Right: Pion mass in 
physical units for the B ensembles together with a constant fit over all data points. 



6. Pressure and Energy Density 

The evaluation of the pressure from the integral technique proceeds by integrating the identity 
^ = r^(^)in temperature along the line of constant physics (LCP): 



T 77-0 t5 



(6.1) 

LCP 



We define the LCP in terms of the pion mass in physical units, which for the smaller mass run is 
shown in the right panel of Fig. |2[ As can be seen it is constant within errors. For the larger mass, 
however, we observe a systematic rise of mps towards larger coupling which amounts to a violation 
of the LCP condition on the level of 10%. 

We perform the integration Eq. ( |0| ) by fitting the available lattice data of ^ to the ansatz 



= exp (-/..f- . \h, + ) , (6-2) 

where t = T /Tq and Tq is a free parameter in the fit. For the fit we use the tree-level con^ected 
data of the trace anomaly that we obtain by normalizing it with the lattice-to-continuum ratio of 



/o{tanh(/i? + /2)} 
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the Stefan-Boltzmann pressure in the free Umit p^^/ psB^ following Ref. iQ]. We check the validity 
of this approach by comparing the continuum limit values as obtained from the corrected as well 
as unconnected data for various temperatures and find compatible results in the majority of cases. 
As can be observed from Figs. ^ and where the thus coiTccted trace anomaly is shown for the 
various available A^^, the correction is efficient and overlays the data from different A^t- 

In order to account for the large errors at small temperatures we perform fits of Eq. ( |6^ ) to the 
upper and lower 1-a deviations and keep the resulting difference as the error of the interpolation. 
For the B (D) ensembles we have fitted data from Nt = 8, 10 and 12 (A^^ = 8 and 10) simultaneously 
and obtain acceptable fits in both cases. We subsequently integrate the interpolation curve numer- 
ically in temperature. The integration constant po in Eq. ( |0| ) has been set to zero in the present 
evaluation. The (yet preliminary) results for the pressure and energy density are shown in the right 
panels of Figs. ^ and ^. 




500 700 
T [MeV] 



interpolation 
3V^ = 10 i 



700 



900 



1 1.5 2 2.5 3 



(e - 3p)/T* 

e/T" C>_- 

Sp/T" c: 

3PSB — 




T [MeV] 



200 300 400 500 600 700 800 

T [MeV] 



Figure 3: Left: The trace anomaly for the D ensembles (see text) obtained for different values of the 



temporal extent Nr- Middle: The trace anomaly after tree-level correction with a fit of Eq. (p.2|). Riglit: 
Preliminary results for the pressure and the energy density. 
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Figure 4: The same as in Fig. |]but for the B ensembles (see text). 



7. Conclusions 

We have presented results for the mass dependence of the pseudo-critical temperature for 
several small values of the pion mass in the range of mps ~ 320 MeV and mps — 470 MeV in a 

'We use Psb/Psb =2.586,1.634,1.265,1.134,1.084 for A^^ = 4,6,8, 10, 12, respectively as obtained in Ref. [0. 
The dependence on the mass (jxT) is found to be very mild and below 1 % such that the same correction factor is used 
for all temperatures. 
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setup with Nf = 2 Wilson twisted mass quarks at = 12. The compaiison of different scenarios 
in the chiral limit is so far inconclusive at the present masses. Further, we have presented first, 
yet preliminary, results of our ongoing project aiming at the determination of the EoS. The trace 
anomaly has been computed for two values of the pseudoscalar mass of about 400 and 700 MeV 
and has been tree-level conected. The pressure has been calculated from the integral method using 
a smooth interpolation formula fitted to the corrected trace anomaly. 
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